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We study ultracold collisions in fermionic ytterbium by precisely measuring the energy shifts 
they impart on the atom's internal clock states. Exploiting Fermi statistics, we uncover p-wave 
collisions, in both weakly and strongly interacting regimes. With the higher density afforded by 
two-dimensional lattice confinement, we demonstrate that strong interactions can lead to a novel 
suppression of this collision shift. In addition to reducing the systematic errors of lattice clocks, this 
work has application to quantum information and quantum simulation with alkaline- earth atoms. 

Ultracold alkaline-earth atoms trapped in an optical field are rich physical systems and attractive candidates for 
quantum information processing [1-4 , quantum simulation of many-body Hamiltonians [5]-[9], and quantum metrol- 
ogy P^QUM] , In each case, interrogating many atoms simultaneously facilitates high measurement precision, but can 
also yield high atomic density and the potential for atom- atom collisions at lattice sites with multiple atoms. For 
quantum information and simulation, these interactions can be a key feature; for quantum metrology, however, they 
present an undesired complication. In either case, these interactions need to be well understood. 

To limit interactions in lattice clocks, the use of ultracold, spin-polarized fermions was proposed to exploit the Fermi 
suppression of 5-wave collisions while freezing out higher partial-wave contributions. However, small collision shifts 
have been measured in fermionic ^^Sr pT [ IT5HI9] and ^^^Yb [12]. Aided by the quantum statistics that govern the 
interactions of these fermionic atoms, we present a complete picture of the cold collisions in the Yb lattice clock by 
performing measurements with state-of-the-art precision together with a quantitative theoretical model. While with 
Sr it was found that 5-wave collisions can occur in the presence of excitation inhomogeneity, with Yb we highlight 
here the important role that p-wave collisions can play in lattice clock systems. Moreover, we demonstrate techniques 
for canceling the collision shift that could be used to vastly reduce the clock uncertainty. 
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FIG. 1: (a) Energy level diagram for Yb. (b) Schematic of the two lattices. Arrows indicate optical polarizations and magnetic 
field directions with respect to gravity. At the left, the Ramsey pulse sequence is entering: two pulses of times ti,2 separated 
by dark time T (not to scale). An inset shows a few 2-D lattice sites with to 2 atoms per site; two atoms in one of the sites 
occupy axial motional states ui and 712. (c) Energy level diagram for two atoms in the rotating frame: three triplet states and 
one singlet state, with interactions V and U, as in Eq. (1). 

To determine the nature of the collisions, here we use two-pulse Ramsey spectroscopy [20] (Fig. 1). Provided 
the Ramsey pulse times ti^2 are short compared to the dark time T, the vast majority of the collisions occur while 
the population is not being simultaneously driven by the laser field, simplifying the analysis of the interactions. 
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The Ramsey technique also ehminates a strong laser-detuning-dependence on the excitation evolution and thus the 
interaction dynamics. Finally, the Ramsey scheme offers the possibility to explore proposals for cancelation of the 
cold collision shift by tailoring the Ramsey pulses [16]. 

In addition to measurements of the collision shift with atoms confined in a one-dimensional (1-D) optical lattice, we 
also show results from a two-dimensional (2-D) lattice, which offers several benefits. First, with strong confinement 
in all but one dimension, the collisions can be treated with a 1-D model. Second, the higher number of lattice sites 
in a 2-D lattice reduces the lattice's overfilling (i.e., many atoms per site). Finally, the 2-D lattice offers stronger 
interactions at any lattice sites that are doubly occupied. 

To model the collisions in a 2-D lattice, we begin by considering a simple case: two atoms in the same lattice site, 
populating axial vibrational modes ni and 722 and the lowest transverse band [4||5l|T7]. Assuming the vibrational 
quantum numbers are conserved during the collisions and laser interrogation, the Hamiltonian for the two-atom system 
can be written in a four-state basis set: \gg), |ee), {\eg) + \ge)) /^/2 ("triplet states") and {\eg) — \ge)) j \f2 ("singlet 
state") ^3[l9j. Here g and e are associated with the lowest ^5*0 and ^Pq electronic levels, respectively, which are 
coupled by the probe laser. The Hamiltonian in the rotating frame can be written in this basis as 



-(5 + K 



^i'^2 72 V2 



^2 '^ni,n2 



(1) 



Here S is the detuning from the atomic transition, l^ni,n2 = (^ni + ^n2)/2 is the average Rabi frequency for the 
two atoms, and Aft = {ftrn — ^712) /2 is the difference in Rabi frequency. The dependence of the Rabi frequency on 
the axial vibrational state is caused by any small projection of the probe beam along the axial direction. The terms 
^nfn2 S^^^' respectively, the s- and p-wave interactions between an a = g^e and a (3 = g^e atom [21]. 
Because the atomic population is prepared in a single nuclear-spin state (m/ = 1/2 or —1/2), quantum statistics 
dictates that only the triplet states, which are invariant under particle exchange, are affected by p-wave interactions, 
while the singlet configuration interacts via 5-wave only. 

For short pulses and large Rabi frequencies, we can ignore interaction effects during the pulses. During the Ramsey 
dark time, the Hamiltonian describing the atom dynamics is diagonal in the singlet-triplet basis, and each state 
acquires just a phase. Consequently, after the second pulse is applied, we recover Ramsey fringes with a frequency 
shift. To compare with the 2-D lattice experiment, the shift must be integrated over the atomic distribution within 
an array of singly and doubly occupied lattice sites. Moreover, instead of choosing a specific set of vibrational modes 
{ni, 722}, we numerically calculate the appropriate thermal average over all possible modes [21 . In the 1-D lattice, each 
site is populated by many atoms, so the two-atom model is not directly applicable. However, in the weakly interacting 
regime, a mean-field picture that approximates the many- atom interactions by a sum of pairwise interactions provides 
a fair description, so we use the two- atom Hamiltonian with temperature-dependent, effective interaction parameters 
to model the multi-atom case. We compared this effective model with a numerically calculated N-body model and 
found qualitative agreement. 

Our experimental procedures are similar to those described in p!2l. After two stages of laser cooling (see Fig. 1(a)), 
atoms are trapped by the horizontal or vertical lattice for 1-D confinement, or by both lattices for 2-D operation. In 
the latter case, we filter away any atoms that are not fully confined at the intersection of the beams by adiabatically 
decreasing the power in one dimension to zero before turning it back up, then doing the same for the other lattice. 
Approximately 2.5 x 10^ atoms are trapped in the 1-D lattice, while ^ as many remain in the 2-D lattice after 
filtering. In the 1-D lattice this corresponds to an estimated density of pi = 3 x 10^^ /cm^ and ~ 20 atoms per site; 
for the 2-D lattice, we estimate that 25 % of the atoms are in doubly-occupied sites, for which the effective density is 
/?2 = 4 X 10^^ /cm^, and fewer than 1 % of the atoms are in sites with more than two atoms. The trap frequencies in 
the lattice were typically 50-75 kHz in the strong direction(s) and 300-500 Hz in the weak direction(s). The lattice 
is tuned to the "magic wavelength" near 759 nm, where the two clock states experience identical trapping potentials. 
We offset the frequencies of the two lattice beams by 2 MHz using acousto-optic modulators (AOMs), preventing any 
line-broadening from the vector Stark shift [T2j [22H24] . 

With the atoms loaded in a lattice, we spin-polarize the sample by optical pumping to one of the spin states 
{rriF = il/2) with 556 nm light; impurity in the spin-polarization is below 1 %. The clock light, pre-stabilized to a 
high-finesse optical cavity [25 to be resonant with the ^5*0 ^ ^Po clock transition, is switched on during the Ramsey 
pulses with an AOM. We interleave high- and low-density clock conditions, each with its own set of integrators to lock 
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the clock laser to the atomic transition and to average over both uif = ±1/2 spin states. The collisional frequency 
shift is found by looking at the difference of the correction signals applied to the AOM divided by the difference in 
atom number, which is typically varied by changing the power to the 399 nm slowing beam that opposes the atomic 
beam. These interleaved measurements have an instability of < 1.5 x 10~^^ / ^/r^ for averaging time r in seconds, 
allowing statistical error bars of ~ 20 mHz in just 2000 s. 
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FIG. 2: (a) Collision shift vs excitation fraction, 1-D lattice. Blue (red) points show experimental measurements in a vertical 
(horizontal) lattice with temperature T ^ 10 /iK and {AQ/Q)j- = 0.2. Dashed black line gives an s-wave-only fit {{U^^)t — 
—27V X 1.5 Hz) from the mean-field model. Solid red line gives a p-wave-only fit with {V^^)r = 10{V^^)r = — 27r x 1.1 Hz. Long- 
dashed green line adds to this a small s-wave component ({U^^)t — — 27r x 0.6). (b) Collision shift vs excitation fraction, 2-D 
lattice. Blue (red) points probe along the vertical (horizontal) lattice. Dashed black line is an s-wave-only fit with a~g ~ —25 ao 



(ao the Bohr radius) ; solid red line is a p- wave-only fit with be 



-74 ao and b^e — O.l^e^. The long-dashed green line adds 



to this a small s-wave interaction 



-25 ao. 



We first considered the collisional shift as a function of excitation fraction (i.e., the fraction in |e) during the Ramsey 
dark time). The excitation fraction was varied by changing the Rabi frequency of the Ramsey pulses. The measured 
shift for atoms in a 1-D lattice is shown in Fig. 2(a) (blue and red points), and for atoms in a 2-D lattice in Fig. 2(b). 
For these data, the Ramsey pulse time is ti = ^2 = 1 ms, and the dark time is T = 80 ms. For the 2-D lattice the black 
dashed and solid red curves give the numerically calculated shift using the s-wave scattering length {a~g) and p-wave 
scattering volumes {h^g and 6^^) as fitting parameters. {V^^ is taken to be zero, consistent with prior measurements 
[26]). For the 1-D lattice, the curves are calculated from the mean-field approximation with the effective interaction 
parameters, which are required to be consistent with those used for the 2-D lattice calculations, varied for fitting. 

Because the atoms are prepared in the triplet state \gg)^ s-wave interactions are allowed only in the presence of 
inhomogeneity (A(^ni,n2 0): which transfers population to the singlet state. The collisions observed in ^'^Sr have 
been attributed to this type of interaction p!5HT9] . By contrast, p-wave interactions are fully allowed, provided there 
is sufficient collision energy to overcome the centrifugal barrier (expected to exceed 30 jiK based on calculated van 
der Waals coefficients [27] )• As shown in the figure, the p-wave interaction provides a much better description of 
the experimental data, as the shift induced by pure 5-wave collisions is generally too small and does not exhibit the 
correct dependence on the excitation fraction. The shifts go through zero near an excitation fraction of 0.51 in the 
1-D lattice and 0.4 in the 2-D lattice. Zero-crossings near 0.5 are readily understood if V^^ dominates: by creating 
equal partial densities of ground and excited atoms, the energy shift on the two clock levels is the same, and the 
net shift is canceled. This effect could allow for vast reduction to the collisional shift in the Yb lattice clock. The 
deviation from a zero-crossing at exactly 0.5 in the 1-D case is consistent with a small ee interaction (6^^ = 0.16^^, 
with beg ^ —74. aQ and ao the Bohr radius). 

We investigated tunneling effects by measuring the shifts for both vertically and horizontally oriented 1-D lattices, 
exploiting gravity-induced suppression of the tunneling rate [28 , but we observed no change in the data (Fig. 2(a)). 
We estimate the tunneling rate, thermally averaged over the lowest nine bands of the 1-D horizontal lattice, to be 
several hertz. In the vertical 1-D lattice and the 2-D lattice, this rate is suppressed by more than a factor of ten due to 
the energy offset between adjacent sites (arising from gravity ^ and the Gaussian beam profile). However, the two 
highest bands of the lattice, which are populated with a few percent of the atoms, can have very high tunneling rates 
(several kilohertz) that are not suppressed by the energy offset. Because of this, the fraction of atoms that can tunnel 
during the spectroscopic time are 10 %, 5 %, and 10 % for the horizontal, vertical, and 2-D lattices, respectively. 
We do not see any appreciable difference between the shifts measured in the horizontal and vertical lattices and thus 
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conclude that tunneling does not play a significant role in the collision shifts. 

We use short Ramsey pulses {ti^2 ^ 1 ms) to avoid interaction effects during the pulse. But, for pulses shorter than 
the mean oscillation period in the trap, 1] > cj^, with i the weakest trap direction, and laser-induced mode-changing 
collisions are not necessarily suppressed. Nevertheless, we ruled out the relevance of those processes by varying the 
pulse duration over a factor of ten without observing any substantial modification to the measured collision shifts [21j . 
We looked for dependence of the collision shifts on the second pulse area [16], but found no significant dependencies. 

To further rule out 5-wave interactions, we misaligned the probe beam to couple more strongly to the weak con- 
finement axis of the lattice trap (Fig. 3(a)). Doing so introduces greater excitation inhomogeneity from the Ramsey 
pulse (in this case, up to a factor of 2.4) because the atoms are not tightly confined along this axis [15 . We expect 
the 5-wave shift to depend quadratically on the inhomogeneity, yet the frequency shifts show no such dependence. 
This insensitivity is well explained by p-wave interactions, which depend only weakly on inhomogeneity at these levels 
(decreasing slightly as more population transfers to the singlet state). A small but non-zero s-wave interaction could 
balance this effect and may help explain the complete lack of dependence, as shown in the theory curves in Fig. 3(a). 
The green long-dashed lines in Fig. 2(a,b) also show that adding a small but non-zero 5- wave interaction is consistent 
with the observed collision shifts. Still, all of these considerations indicate that p-wave interactions play the dominant 
role in the cold collisions of ^^^Yb. 
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FIG. 3: (a) Collision shift vs probe misalignment angle (vertical 1-D lattice) for constant excitation fraction 0.12. Using the 
same parameters as Fig. 2(a), the dashed black line gives an s- wave-only fit, solid red line gives a p- wave-only fit, and the 
long-dashed green line has s- and p-wave terms. In the well-aligned case (0 mrad) there is a residual effective misalignment 
of ~ 5 mrad due to the imperfect overlap between lattice and probe beams, (b) Collision shift vs Ramsey time, 1-D lattice, 
for excitation fraction 0.18. Using the same parameters as Fig. 2(a), the solid-line gives a fit from the mean-field model, 
(c) Collision shift vs Ramsey time, 2-D lattice, for excitation fraction 0.19 =b .03. The shift crosses zero due to the periodic 
dependence of the shift on collisional phase, and is a signature of strong interactions. The model calculations (shaded region) 
use the same parameters as Fig. 2(b) for an excitation fraction range 0.19 di .03. (d) Asymmetric Rabi spectrum, 2-D lattice. 
The solid line is the prediction of the model, using the number of doubly occupied sites as a fitting parameter. 

Strong interactions emerge in the regime V"^T > 1. A key observation revealing the operation of the 2-D lattice 
clock in the regime V^^T ^ V^^T > 1 is the zero-crossing of the collision shift at a lower excitation fraction of 
0.4, which deviates from the crossing at 0.5 predicted by the weakly interacting expression of the shift [21]. The 
interaction strength also introduces additional dependencies on the Ramsey dark time. With weak interactions, the 
collision shift is independent of T, but with strong interactions the clock shift decays with increasing T due to the 
shift's sinusoidal dependence on scattering phase [2TJ|29]. We investigated this experimentally by varying the dark 
time T and measuring collision shifts in the 1-D lattice (Fig. 3(b)), where the shift scales weakly with T, and the 
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2-D lattice (Fig. 3(c)), where the shift is strongly damped towards zero with increasing T. Yet a third signature of 
strong interactions is significant asymmetry in the clock transition spectrum. In Fig. 3(d) we show a Rabi spectrum 
(t = 120 ms), taken under high density operation in the 2-D lattice, which shows an additional feature on the red 
side {S < 0) of resonance. This asymmetry is density-dependent and barely observable in the 1-D lattice. In the 2-D 
lattice, the interactions are sufficiently strong {V^^t > 1) to introduce these asymmetric lineshape features beyond 
the transition linewidth. With yet higher density, it may be possible to spectrally resolve three features, one each 
for the 5-wave-interacting singlet, the p- wave-interacting triplets, and the non-interacting atoms in singly occupied 
lattice sites. Interaction-induced sidebands were recently reported in ^ and may be useful for quantum simulation 
applications. 

In this Letter we have shown evidence for p-wave interactions in ultracold Yb confined in an optical lattice. Although 
lower atomic temperature yields reduced tunneling through the p-wave barrier, and thus a lower scattering cross- 
section, it also increases the atomic density of the confined atoms. For this reason, both s- and p-wave interactions 
may be potentially relevant for all optical lattice clock systems. Using the dependence of the measured shift on 
excitation fraction and Ramsey dark time, we have observed zero-crossings in the measured frequency shifts, which 
provide the metrological means to reduce the shifts to nearly negligible levels. 
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Supplementary Material 



TWO ATOM MODEL 



Here we consider two nuclear-spin-polarized fermionic atoms interacting via s-wave and p-wave channels. The 
atoms are confined in a tube with frozen transverse degrees of freedom (only the lowest vibrational transverse mode is 
populated). Along the tube direction, Z, there is a weak harmonic confinement with frequency ujz^ and we will first 
assume the two atoms populate the vibrational modes ni and 712- In this case there are only four states spanning the 
Hilbert space: the triplet states \gg)^ |ee), {\eg) + \ge))/\/2^ and the singlet, {\eg) — \ge))/^/2. Here the convention 
used is that the left atom populates mode ni and the right atom mode n2. 

In the presence of a laser field with wave vector k = kyY -\- kzZ and detuned from the atom transition 
frequency by ^, the two-atom Hamiltonian in the rotating frame reduces to Eq. 1 in the main text. Here 

Q.n = ^o^n(^|)^o(^Y)e~*^^^^^^^/^, with rji = ^^uj- Lamb-Dicke parameter along the i direction, Vt^ the 
bare Rabi frequency, and the Laguerre polynomial, m is the atom mass. The 5-wave and p-wave interaction 

parameters are V^^.^^., = 4^^^^^^^a- 5'^^.,^ and V^f^^,, = 12^^^^^^^6^^^Pn,-,n^., • Here 5^^.,^., and Pn,.,n^., are 
geometric terms that take into account the spatial overlap between the atomic wavefunctions of colliding atoms in 

I ci±±vj. ± Yin' — / ; 5 

where Hn are Hermite polynomials. The dependence of V^^ and U^^ on the vibrational mode encapsulates the 
temperature dependence of the interactions. 

We ignore interactions during the pulses, and consequently the number of excited atoms after the first pulse is 

During the dark time, the Hamiltonian is diagonal and each state acquires just a phase. After the second pulse, the 
excited state population is 

N^,,n,{tiM = An,,n, ^^n„n, cos[{S - 27rA<^,jT], (3) 

with An^^ri2 an overall offset, Mm,n2 > the fringe amplitude, and Ai/^^ ^^le frequency shift. These quantities can 
be computed analytically. 

In the 2-D lattice system an array of isolated tubes is populated with mainly one and two atoms per tube. If TVq 
tubes are singly occupied and Ni tubes doubly occupied, then 



modes rij and n./. They are given by Snn' = ^ i — ^^'-^-'^^'^^ 1^ and Pnn' = 



with 



(ii, i2) = + cos[,5T] {A\^^^^ + + sin[,5T]S^^,„,, (4) 



^1 . . . . 

A^^ = [1 ~ (^^l) (^^2) '^i" (^1) 'Si" (^2) (^2T (y^^ - U^i^y^ - cos (a^^) cos (A6I2) cos {e\^ cos (5) 

1 ^0 

+ - XI ~ (^^1) (^^2) (^1) sin (^2) - cos (a0^) cos (a0|) cos {e\^ cos (^2)]! 

2 i=l 

1 ^1 
^ i = l 

- cos (TC^) \cQS (A6»2) sin sin {d\^ sin (TB^) + sin (^2A0l^ sin (A6'2) cos (^6*2) sin (TD^)^ 

+ 2 sin (A^f) sin (Ae*) cos (e^) cos (e^) sin (TC^) cos (TI)^) + cos (Ae^) cos (a^) sin (^f) sin (0^) [sin (t (vf^ - V^^)) + (^^ (v^^ _ l^SS))] ; 



1 

2 , 

1 'JU 



(7) 



<,n2 = -E («) 



cos (A6I2) sin (^el^ sin (^82) [2 cos (^6*^) sin (TB^) sin (TC^) + cos (Ae*^) (^cos (^T (^vf ® - V-^^)) + cos (^T (^V®^ - V^^^))] + 
sin (a4) cos (4) [sin (2A0f ) sin (TC^) sin (TI?^) + sin (a^^) cos (^f ) (cos (t (^U^^ - V/^)) + cos (t ([/^^S - V^S)))] . 

In the above equations, the dependence of the various parameters on the modes {ni, 712} is omitted for simplicity but 
implied. The index i runs over the set of doubly- and singly-occupied tubes. {Vf^ — V^^) = 2Q, V^^ + V^^^) = 
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2Bi, and {—2U^^ + Vf^ + F-^^) = 21^^. Due to the Gaussian profile of the laser beams, the trapping confinement varies 
from tube to tube. This variation gives rise to tube-dependent interaction parameters as well as tube-dependent Rabi 
frequencies $1 = tg^^ and AOl = tsAQ\ Here 5 = 1,2. 

If the atomic population is initially prepared in the excited state (|e)) instead of the ground state (|^)), then 

N^,,nAh,t2) = ^i^„„.+cos[5T]«_„^+^2g,„J+sin[5T]S^^^,„^, (9) 

with 

^^1 ,r.2 = + (^^i) (^4) (^l) (^2) (V/^S - U^^)) + cos (a^^) cos (a^^) cos (a^) cos (o^)] (10) 

1 ^0 

+ - E [1 + (a^I) sin (a^^) sin (ol) sin (o^) + cos (a^^) cos (a^^) cos (ol) cos (e^)]; 
^ i = l 

1 ^1 
^ i = l 

cos (a4) sin [sin (2^^) sin (TB,) cos (TC,) + cos (a^^) sin (^^) (sin (t (v^ - V^^)) + sin (t (v/^^ - vf^)))] + 

sin (A6»2) cos (6*2) [2 sin (a^^) cos (^^^ sin (TC^) cos (TD^) + sin (2A0^) cos (TC^) sin (TD^)^ ; 

^ i = l 

cos (A6>|) sin (6>f ) sin (6»|) [cos (aS^) (cos (t (vf ^ - V^.^^)) + cos (t {v^^ - ^))) - 2 cos (e^) sin (TB^) sin (TC^)] 

+ sin (a^^) cos (el) [sin (A0^) cos (e^^) (cos (t (c/^^^ - ^)) + cos (t (t/^^S - vf^))) - sin (2Ae^) sin (TC^) sin (TD^)] . 



In both cases the shift is given by 

arctan ( ,_^2^^"Xa^^^ w ) + ""P 
{Av)r = "'-^^^^"""^^"^ , (13) 

with Ot" denoting a thermal average, p is an integer that needs to be chosen so that during the g ^ e {e ^ g) 
interrogation the total number of atoms driven to e {g) reaches a maximum value at 5 = 27t{Au)j-^ instead of a 
minimum. It also ensures that the shift is a smooth function of T. If p is not correctly chosen the shift can jump 
discontinuously instead of becoming smoothly displaced outside the first Ramsey fringe. For the experimental regimes 
described here, no discontinuity occurs and we can set p = 0. 

An important point to emphasize is the different dependence of the shift on T in the weakly and strongly interacting 
regimes. In the weakly interacting regime, TV"'^ <C 1 and TU^^ <C 1, the B term provides the leading contribution, 
because it exhibits a linear dependence on interactions B oc (TV"'^, This results in a T-independent collision 

shift. On the other hand, in the strongly interacting regime, (TV"'^, ^ 1, both terms B and A exhibit a 
nontrivial sinusoidal dependence on the various interaction parameters. This implies that as T increases, both B and 



A exhibit faster but bounded oscillations and therefore, according to Eq. 13, on average the shift decays as 1/T. 



INTERACTION SIDEBANDS: RABI SPECTROSCOPY 



We consider now the case in which the system is interrogated during a long single Rabi pulse of duration t. In 
this situation, both interaction- and laser-driven terms must be accounted for at the same time during the dynamical 
evolution. For s- wave-dominated collisions, reaching the strongly interacting regime U^^t > 1 can lead to a suppression 
of the clock shift [171 |T9l . This is not necessarily the case for p- wave-dominated collisions, in which atoms interact 
even in the initially populated triplet manifold. 

Moreover, if {V^^ — V^^)t > 1 and the atoms start in the triplet \gg) state, there is a large gap that they must 
overcome if (5t <C 1 to populate the + \eg)) state. This means that only when 6 — V^^ — V^^ do the two states 

become resonant, and population will be transferred between \gg) and -^{\ge) + \eg)). Since \gg) is coupled to \ee) 

by a second-order process, population of \ee) is energetically suppressed. Only if S = is there a two-photon 

resonance, but even in this case the overall amplitude is small since it is proportional to Cl'^/Bn^^n^. 
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From these considerations, one can write the hneshape of the 2D lattice array as 



ni ,712 



ni ,712 



/ Ni Ni 

\i=l i=l 

Ni No \ 

+ ^ p[V2AQl,^^^^ , 5 - Ul\%^ + J + ^ , J] ) (14) 



with y\ = ^2j^y2 sin^[yx^-r^t/2]. For comparison with experiment, we need to average over the tube array and 
evaluate thermal averages. If it were possible to resolve the various peaks we would expect the line shape to have three 
single-photon resonances: One at ^ = coming from single occupied tubes, one p-wave related at (5 ex (V^^ — V^^)r^ 
and one at ^ ex {U^^ — V^^)r due to 5- wave collisions. 

To obtain the excitation fraction and the interaction sidebands, we need to properly integrate over the spatial 
degrees of freedom of the lattice. The characteristic harmonic oscillator frequencies uoi {i = X, F, Z) of each lattice 
site vary smoothly with the position of the lattice site modifying the Lamb-Dicke and interaction parameters. The 
variation of the trapping frequencies is obtained from the lattice potential 

V{X, r, Z) = -Vx,o cos^ f ^ ) e- w2 _ Vy^o cos^ f ^ ) e' , (15) 

where a is the lattice spacing, W ~35 jam is the beam waist, and Vx^o and VV,o are the lattice depths. We extract 
the characteristic harmonic oscillator frequencies of a given site by expanding the lattice potential at the center of 
it. Assuming a uniform occupation of the sites in the center of the lattice, we transform the spatial integration into 
an integration over the relevant range of trapping frequencies weighted by the number of sites with such trapping 
frequencies. For each set of trapping frequencies, a thermal average is carried out over both the harmonic oscillator 
modes and lattice vibrational modes. 

This procedure is significantly simplified by noting that the main correction of the spatial average comes from the 
changes in the interaction terms. Thus, changes in the Rabi frequencies are accounted for only by using an average 
value, allowing us to transform the integration over trapping frequencies into an integration over interaction strengths. 



OTHER CONSIDERATIONS 



We have restricted so far our analysis to an effective one-dimensional system with at most two atoms per tube. 
For more than two atoms per tube the spin model is no longer diagonal in the collective spin basis, since the spin 
Hamiltonian is not SU(2) symmetric. This implies that for more than two atoms and outside the weakly interacting 
regime our calculations are not exact. However, we have performed numerical calculations for up to five atoms per tube 
that show that the two-atom model can qualitatively predict the many-atom case, but with renormalized interaction 
parameters. 

In the 1-D lattice there are two weakly confined directions and multiple atoms per lattice site. Under these 
conditions, interaction-induced mode-changing collisions, which are not accounted for in our two-atom model [17], 
have to be included. We have studied the role of those processes by numerical evaluation of a more general multi-mode 
Hamiltonian and found that it mainly leads to a renormalization of the model parameters. All of these considerations 
justify the validity of the two-atom model for the description of the 1-D lattice. 



EXPERIMENTAL DATA ON VARIOUS PULSE DURATIONS 



We use short Ramsey pulses (ti^2 ^ 1ms) to avoid interaction effects during the pulse. However, for pulses shorter 
than the mean oscillation period in the trap, Q > Ui with i the weakest trap direction, and laser-induced mode- 
changing collisions (not included in the two-atom model) are not necessarily suppressed. To explore the role of 
laser-induced mode-changing collisions, we measured the collision shift as function of excitation fraction in the 2-D 
lattice for various pulse durations. As shown in Figure [4j we varied the pulse time ti^2 by one order of magnitude but 
found no substantial change in the measured shifts. Thus we conclude that these processes do not play an important 
role. 
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FIG. 4: We study laser-driven mode-changing collisions by measuring the collisional frequency shift vs. excitation fraction for 
various Ramsey pulse times in the 2-D lattice. The blue, red, and green points correspond to ti,2 = {0.5, 1, 5} ms, respectively. 
Because we observe no substantial trend, we conclude that these processes do not play an important role. 



